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ABSTRACT 

A mathematical model is proposed for closing or 
matnematical ly completing the system of equations 
wnicn descriDes the time-average flew field through 
the Dlaae passages of multistage turoomachinery . 

Tnese equations referred to as tne average-passage 
equation system govern a conceptual model which has 
proven useful in turoomachinery aerodynamic design 
and analysis. Tne closure model is developed so as to 
insure a consistency between tnese equations and the 
axisymmetric tnrough-flow equations. The closure 
model was incorporated into a computer core for use in 
simulating tne flow field about a hign-speed counter- 
rotating propeller and a high-speed fan stage. 

Results from tnese simulations are presented. 

INTRODUCTION 

Engineers nave <ong recognized tne difficulty 
associatea with adopting a “First Principle" approach 
Dased on directly solving the Navier-Stokes equations 
for the purpose of designing (or analyzing) vehicles 
wnicn operate in high Reynolds number turbulent flows. 
However, numerous examples exist, suen as turbomachin- 
ery olading, aircraft wings and bodies, inlets and 
nozzles, wnich clearly show that models which describe 
an “averaged" flow state can be used to design aero- 
dynamic venicles and provide answers to many aero- 
dynamic proolems. In Doth external and internal 
aerodynamics, the “averaged" state most often modeled 
is one in wnich tne flow appears steady, in general, 
tne numoer of equations associated with this 
"av'.ageo" flow representation does not equal the 
numoer of unknowns. Tne problem of mathematically 
completing this system of equations so that they may 
oe solved is referred to as tne closure proolem. The 
flow moaels associated witn the completed system of 
equations must De considered semi-empirical for they 
rely neavily on empirical correlations to introduce 
tne effects of turoulent motion ano, in the case of 
turbomachinery, the additional effects of unsteadiness 
and spatial nonunifor nities into tnese “averaged" flow 
representations. For nonturbomachinery application, 
tne equation governing such a flow is the familiar 


Reynolds-averaged Navier-Stokes equation. In general, 
tne lengtn scales associated with this equation are 
sufficiently restricted so as to make tnem amenable to 
nu:nerical simulation. Indeed tnere is considerable 
accivitv tnese days in the external aerodynamic com- 
munity co develop numerical simulators Dased on these 
equations for flows over an entire aircraft. 

For turDomacninery involving more than one blade 
row, the Reynolds-averaged form of tne Navier-Stokes 
equations do not describe a flow which is steady in 
time. On tne contrary, they descriDe a flow whicn is 
nignly unsteady in wnicn Dlade rows are moving rela- 
tive to one anotner, generating disturbances whose 
time scales range from a fraction of wneel-speed to 
many times that of Dlade passing frequency and wnose 
length scales range from tne circumference of the 
macnine to the thickness of the laminar sublayer 
region of the turDulent boundary layers. Simulation 
based on tne Reynolds-averaged Navier-Stokes equations 
are well beyond the capabilities of today's computers 
for all Dut the simplest of multistage geometries. 

Tney also do not govern the conceptual flow model tra- 
ditionally used to design multistage turbomachinery. 
as noted above, mult tage designs are based on flow 
models in wnich the fkw appears steady within each 
Dlade row. In addition, •.lith respect to a given blade 
row, tnese models assume the flow to be spatially 
periodic from one blade passage to anotner. In 
Ref. 1, a mathematical derivation of the equations 
governing this flow was presented. Tnese equations 
were referred to as the average-passage equation 
system. This derivation was carried out for arbitrary 
configurations and clearly showed the relationsnip 
Detween tne Navier-Stokes equations, tneir Reynolds- 
averageo form, and tneir average-passage form. Tne 
closure problem associated with the average-passage 
form of tne Navier-Stokes equations was also identi- 
fied. Tnis work put tne average-passage model on a 
sound matnematical base equivalent to that of the 
Reynolds-averaged Navier-Stokes model. A brief sum- 
mary of that work is presented in tne next section. 

Tne purpose of tne present work is to elaborate 
furtner on the issue of closure for the average- 
passage equation system and to propose a closure model 
for the inviscid form of this equation system. This 
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closure model was used to obtain the results presented 
in Ref. 2. Tnose results and the ones to be presented 
in this work show that the present model appears 
applicable to configurations in wnich the average- 
passage flow field is nearly irrotational between 
blade rows. 

MODEL EQUATION HIERARCHY 

In Fig. 1, a hierarchy of equations are shown 
wnicn can be used to analyze turbomachinery flows. 

Tne Navier-Stokes equations appear at the upper left- 
nano corner of this figure. These equations are 
assumed to provide a complete description of the flow 
field, including a complete description of turbulent 
motions. To use these equations as a basis for simu- 
lating turbomachinery flows requires sufficient com- 
puter capacity to resolve all of the time and length 
scales associated with high Reynolds number flows. In 
addition, since turbomachinery flows are statistically 
nonstationary, a sufficient number of computations 
would nave to be performed over a range of randomly 
chosen initial conditions to insure a stastistical 
steady-state description of the flow. Such simula- 
tions are clearly Deyond the capacity of today's most 
advanced computers. The next box (i.e.. Fig. 1) con- 
tains tne Reynolds-averaged form of the Navier-Stokes 
equations. Tney are derived by ensemble averaging the 
Navier-Stokes equations and hence govern a determin- 
istic description of tne flow field. An illustration 
of tnis description for a two-stage configuration in 
wnicn tne first and second rotors have five and four 
Dlades respectively while the first and second stage 
stators nave four and five blades is presented in 
Fig. 2. The rotors rotate relative to the stators, 
and, therefore, the flow will be unsteady in either 
tne rotor or stator frame of reference. As noted 
previously, the time scales associated with this 
unsteady flow are quite diverse, which makes simula- 
tion for all out tne most simple of geometries beyond 
the capabilities of today's computers. The closure 
problem associated with these equations requires the 
modeling of tne familiar Reynolds stress and energy 
correlations. It is by means of these correlations 
that the "average" effects of random fluctuations in 
momentum and energy of a fluid particle are introduced 
into the equations governing the deterministic flow 
field. 

Tne third box from the left in Fig. 1 represents 
the time-averaged form of the Reynolds-averaged 
Navier-Stoxes equations. These equations govern the 
time-averaged flow field as viewed by an observer 
whose frame of reference is fixed to a given blade 
row. An illustration of this description for the 
two-stage configuration used to illustrate tne 
Reynolds-averaged flow model is also presented in 
Fig. 2. All rotating blade rows have a unique time- 
averagea flow field associated with tnem. In a 
similar fashion, all nonrotating blade rows have tneir 
own time-averaged flow field representation. These 
two flow fields are not the same. For both flow 
fields, the blade rows which rotate relative to the 
olade rows wnicn are stationary (i.e., witn respect 
to one anotner) appear smeared. Their physical 
appearance is very similar to what one observes wnen 
viewing a high-speed propeller. Within the context of 
the time-averaged flow description, these smeared 
Dlade rows are replaced by actuator ducts (i.e., 
actuator ducts of finite thickness). These ducts are 
represented oy a oody-force distribution which can add 
or extract energy from the flow. In addition, the 
time-average flow equations contain correlations 


between time-varying flow variables. These correla- 
tions arise because the Reynolds-averaged Navier- 
Stokes equation is nonlinear. These correlations 
represent the time-average of tne fluctuating density 
field and products of the fluctuating velocity field 
as well as the time-average of the fluctuating 
density, fluctuating velocity, and fluctuating total 
enthalpy field, it is through these correlations 
that the "averaged" effect of the relevant unsteady 
physical phenomena is introduced into the time- 
averaged representation. The modeling of the body 
forces and energy sources associated with the smeared 
blade rows and the temporal correlations, plus the 
modeling of the time-averaged Reynolds stresses, forms 
tne closure problem associated with the time-averaged 
equations. Finally, it should be noted that, for a 
single-stage configuration, the time-averaged flow 
field associated with either olade row will be spa- 
tially periodic over the pitcn of that blade row. 

Thus, if the closure issue associated with the time- 
averaged representation can be addressed without over- 
due complexity, it should be feasible to conduct a 
simulation based on this flow model for a single 
stage. 

For a multistage configuration in which the 
number of rotor blades differ from rotor to rotor, or 
for wnich the numoer stator blades differs from 

stator to stator, the me-average flow field will 

not, in general, be spatially periodic over the pitch 
of any given blade row. An averaging-procedure may 
oe introduced which transforms this spatially aperi- 
odic flow field into one that is periodic over the 
pitch of a given blade row. Tne resulting flow field 
is referred to as tne average-passage flow and appears 
in the fourth box from the left in Fig. 1. Each blade 
row in a multistage machine has associated with it an 
average-passage flow field. An illustration of this 
description is shown in Fig. 2. For the two-stage 
machine under consideration there exists four average- 
passage flow descriptions due to the number of blades 
assigned to each wheel. The geometry of neighboring 
blade rows (rotating and stationary) for which the 
blade count is not an integral multiple of the blade 
row of interest, and are stationary relative to this 
Dlade row, appear smeared in this flow description. 
Their appearance is similar to that of the rotating 
blade rows in the time-averaged flow description. It 
snould De noted that all of the blade rows which 
rotate relative to the blade row of interest appear 
smeared, since the average-passage description is also 
a time-averaged description. The four average-passage 
flows illustrated in Fig. 2 are coupled to one another 
through a system of body forces, energy sources, and 
temporal and spatial correlations. Tne closure prob- 
lem associated with this flow description consists of 
developing matnematical expressions for the spatial 
and temporal correlations in addition to the body 
forces and energy sources. These correlations intro- 
duce tne transport on the “average" of momentum and 
energy between the time-averaged representation and 
the average-passage representation. 

Many analyses currently used to analyze multi- 
blade row turbomachinery involve iterating between a 
meridional flow analysis and a blade-to-blade anal- 
ysis. Within the context of the present discussion, 
these analyses may be viewed in one of two ways. They 
may be thought of as attempting to describe the 
average-passage flow field. If one gives these anal- 
yses this interpretation, then one immediately notes 
that their derivation lacks matnematical rigor. As a 
result, the closure problem associated with the 
average-passage representation is never addressed. 
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for it is completely overlooked. On the other hand, 
one may Interpret these as axisymmetric analyses In 
wnicn the blade-to-blade solution, along with some 
empirical correlations, is used to close the merid- 
ional flow equations. In this case, these analyses 
are rigorous because the closure problem that being 
tne closure of the axisymmetric representation is 
generally clearly defined. However one wishes to 
interpret these analyses, one must be impressed with 
the degree of accuracy with which they predict the 
axisymmetric flow field in the neighborhood of design 
conditions. As one moves away from the neighborhood 
of the design point, however, the validity of these 
analyses appears to degenerate quickly. This dis- 
agreement is thought to be due to the inability of the 
blade-to-blade model to properly account for large 
spanwise migration of flow which occurs at these off- 
design conditions. To analyze such situations, a true 
three-dimensional analysis is needed. The average- 
passage model provides a framework for developing such 
an analysis, as illustrated by the work presented in 
Ref. 2. The accuracy of such simulation will, of 
course, depend upon the validity of the closure model 
used in the simulation. 

Tne next box in Fig. 1 represents the axisym- 
metric flow model, which is tne mainstay of many turbo- 
macninery design systems. Tne field equations for 
this model can oe derived oy tangentially averaging 
the average-passage equation system. An illustra- 
tion of the geometry associated witn this representa- 
tion is also provided in Fig. 2. Each average-passage 
flow model can be related to an axisymmetric model. 

Tne equations governing these four axisymmetric models 
must be equal to one another, for there can only be 
one axisymmetric or througn-flow representation of the 
flow field witnin a multiblade row configuration. 

Tne average-passage equations thus define the three- 
dimensional passage flows naving a common axisymmetric 
flow description. All of the blade rows within the 
axisymmetric description appear smeared and are mathe- 
matically replaced by actuator ducts. These ducts 
exert a force on tne fluid which may add or extract 
energy from tne flow. There may also be energy 
sources or sinks witnin tne ducts wnich are associated 
witn blade heat transfer. Over tne years, numerous 
publications have appeared which dealt with modeling 
tnese forces and the energy sources. Quite often they 
are estimated from cascade or blade-to-blade analyses 
tempered oy empirical correlations. In addition, the 
axisymmetric or through-flow equations contain corre- 
lations between temporal varying flow variables as 
well as correlations between spatial varying flow 
variaoles. Tnese correlations introduce on the 
“average" the effect of radial transport of momentum 
and energy from the average-passage representation. 

Only very recently nave models for these correlations 
appeared in the open literature. Sehra (3) was one of 
the first to attempt to incorporate tnese’correlations 
into a througn-flow code. His correlation model was 
Dased on data obtained from a high-speed isolated 
rotor test. Jennions (4) modeled tnese correlations 
using results from an inviscid olade-to-olade anal- 
ysis. He was able to develop an iterative procedure 
for incorporating these correlations into a tnrough- 
flow analysis. Finally, the Adkins and Smith (5) 
model for accounting for the effects of the spanwise 
mixing in multistage machinery, may be thought of as 
an attempt at modeling the correlations which appear 
in the axisymmetric model. Tne last box in Fig. 1 
represents a quasi-one-dimensional equation system. 
Tnese equations result from averaging the axi symmetric 
equation over the span of the flow annulus. Tnis 
equation system is often used in engine stability 


studies and in preliminary design to establish the 
flow properties along the pitch-line of a machine. 
Closure of this system of equations can be quite 
involved. It requires models for the blade forces, 
energy sources, spatial and temporal correlations 
associated witn the o 1 ade-to-b 1 ade flow field, as 
well as a model for the force exerted by the casing 
on tne flow. 

The flow models identifieo in Fig. (1) are by no 
means complete, nor were they ever intended to be 
complete. The purpose of this figure was to illus- 
trate symbolically the connection between a hierarchy 
of equations associated with turbomachinery aero- 
dynamics. It is hoped that the rational derivation 
of tne average-passage equation system will ultimately 
'ead to the development of three-dimensional viscous 
computer codes for multistage configurations. Such 
codes will enhance our ability to analyze turbo- 
inachinery flows, especially at off-design conditions. 
For it is our inability to accurately predict off- 
oesign performance of multistage machinery which is 
often tne major contributor to their high development 
costs and not proDlems associated with poor design 
performance. In the next section, the closure model 
associated witn tne inviscid form of the average- 
passage equation system will be developed. 

THE CLOSURE PROBLEM 

For simplicity we shall only address the closure 
proDlem associated with solving the inviscid form of 
tne average-passage equation system as it pertains to 
a single stage. A solution to the corresponding 
multistage problem can be obtained by a direct exten- 
sion of the analysis wnicn follows. For a single- 
stage configuration, each blade row has associated 
witn it an average-passage equation system. As noted 
in the previous section, the dependence of the flow 
througn the first blade row upon that through the 
second is introduced by means of a body force, energy 
source, and time-average correlations between fluctu- 
ating flow variables. Likewise, a corresponding 
dependency exists between the flow through the second 
olade row and that tnrough the first. In Ref. 1, the 
body force and energy source which appear in the 
inviscid form of tne average-passage equation system 
were snown to depend upon tne ensemble-averaged pres- 
sure. This "averaged" pressure was estimated from 
samples of tne pressure field taken over a period of 
one revolution of tne wneel recorded at tne instant a 
olade passes an observer whose frame of reference is 
fixed to that of tne blade row of interest. If one 
assumes tne average-passage flows of the two blade 
rows to oe nearly irrotational outside of tne blade 
passage region, then tnis ensemole-averaged pressure 
is nearly equal to the average-passage pressure dis- 
tribution on the surface of the neignboring blade row. 
Hence, from tne solution for the first blade row, one 
can estimate the body force and energy source which 
appear in the equations for the second blade row. In 
a similar fashion, one may estimate the body force and 
energy source which appear in the equations for tne 
first blade row from a solution to the corresponding 
equations for the second blade row. 

The remaining terms which muse be estimated are 
tne temporal correlations associated with the time- 
varying flow field. The origin and nature of these 
correlations were discussed in the previous section. 

To develop a model for these correlations, we 
decompose the absolute velocity field, it, according 
to the equation 
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Y (r,e,z,t) ."v^hr.i) 

♦ (r.e-Qjt.z) - v^ AX * (r,z)j 

♦ |V 2) (r,e-n 2 t,z) r,z)j 

♦ v (3) (r,e,z,t) (1) 

where tK ax ) represents the axisymmetric velocity 
component, the time-averaged absolute velocity 

field as observed in a frame of reference fixed to the 
first Dlade row, the corresponding velocity 

field observed in a frame of reference fixed to the 

second blade row, and the component of 

velocity which is unsteady in either frame of refer- 
ence. Tne remaining variaDles which appear in Eq. (1) 
are the cylindrical coordinates r,u,z, time t, and 
tne rotational speed of the first and second blade 
rows az- In a similar fashion, the total 
entnalpy, H, measured in the absolute frame of refer- 
ence can be decomposed according to the tqua.ion 

H (r.e.z.t) - H^ AX )(r,z) 

+ [h ( 1) (r.e-s^t.z) - H (AX) (r,z) j 

+ (r,e-n 2 t,z) - H^ AX * (r,z)J 

♦ H (3) (r.e.z.t) (2) 

Fixed in the frame of reference of the first blade row, 
tne velocity field will appear steady in time, 

wnile the components "v' 2 ' and v”) w ni appear to be 
unsteady. If we define the velocity component 
v"(r,e,z,t) as 


■$“ (r.e.z.t) i v (2) (r,e-a 2 t,z) -v^ AX ^ (r.z) (3) 

tne correlations whicn appear in tne average-passage 
momentum equations associated with the first blade row 
are obtained by forming the time-average of the 
product of the fluid density and (?" ♦ ?'•*'), (1) . 

Tne result is 
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(4) 


where the subscripts i,j take on the values of 1, 2, 
and 3. On the right-hand of Eq. (4), these subscripts 
are used to denote the axial, tangential, and radial 
velocity components respectively. The variable p is 
the fluid density, and the over-bar represents the 
time-average of the variables which appear beneath it. 
Thus for i . 1 and j - 2, denotes the temporal 
correlation between the density p, and the product of 
tne axial and tangential components of the fluctuating 
velocity field. For low Mach number flows in which 
tne den sity may be assumed constant, the correlation 
p Vi vj will be independent of tangential position 
e, since v? is spatially periodic over the pitch of 
tne second blade row. This correlation is thus asso- 
ciated with tne transport on the “average" of momentum 
across the axisymmetric stream surfaces. Tne remain- 
ing correlations whicn appear in Eq. (4), however. 


will be functions of e if is spatially 

aperiodic over the pitch of the second blade row. In 
general this will be the case. 

Based on the arguments used to derive Eq. (4) and 
the analysis presented in Ref. (1), the correlations 
which appear in tne energy equation are obtained by 
forming the time-average of the product of p, 

H“ + Hp) and v" + v”'. The result is 

Q i - pH" v^ ♦ pH" v| 5) + pH^ v^ + ^IT 


(5) 


wnere 


H" (r.e.z.t) = (r,e-n 2 t,z) - (r.z) 


(b) 

The first correlation in Eq. (b) is independent 
-f e if the fluid density is constant. This is 
the result of H" and vi“ being spatially peri- 
odic over the pitch of tne second blade row. In 
general, the remaining correlations will not exhibit 
this behavior. As a result, the total enthalpy asso- 
ciated with the average-passage flow field for a 
multiolade row configuration wi’l be nonuniform in 
the tangential direction. Ken ebrock and Mikolajczak 
(b) were the first to attempt to analyze the fluid 
mecnanics associated with this phenomena. They attri- 
buted it to the transport of excess total temperature 
of a fluid particle in a rotor wake across the stator 
passage. Their analysis of this process was based on 
kinematics. Although the present work makes no 
attempt at developing an alternative model of this 
phenomena, it does suggest that it is associated with 
the dynamics of stator-blade rotor-wake interaction. 

For an inviscid nearly irrotational flow, the 
magnitude of tne unsteady component w jn &e 

comparable or less than the magnitude of "v", except 
for regions near blade leading edges. In particular, 
in regions where the body force and tne energy source 
are finite, the correlation associated with v" 

(i.e., Eqs. (4) and (5)) will be significantly larger 
tnan those associated witn 'v^'. For this reason 
we assume tnat the correlations in Eqs. (4) and (5) 
associated witn the unsteady velocity component Y3) 
can be neglected. As a result, the correlation Rjj 
ana the correlation Qj can be directly evaluated 
from the average-passage solutions. For a stage, this 
implies that the flow field tnrougn both blade rows 
must be evaluated simultaneously. 

To incorporate the suggested closure model into 
a numerical simulation, one may envision a two-tier 
iteration procedure as depicted in Fig. 3. In the 
inner loop, the body forces, energy sources and cor- 
relations are frozen. An average-passage flow field 
is evaluated based on tne value of these quantities 
and the imposed boundary conditions. In tne outer 
loop, the body forces, energy sources, and correla- 
tions are updated based on the converged inner loop 
solutions. Me must update these terms, as previously 
noted, in a manner which yields a unique axisymmetric 
representation of tne flow field through the machine. 
Tnis will insure that the average-passage representa- 
tion of the flow is consistent witn the axisymmetric 
representation. 

Tne equations to be solved in the outer loop may 
be derived starting from tne equations of motion 
expressed in the vector form (2). 
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L(U„) ♦ y K(U n )dV ♦ y* S n _! dV - 0 (7) 

The operator L(U n ) in tnis equation represents 
tne net flux of mass (i.e., continuity of flow through 
a control volume), axial and radial momentum, angular 
momentum, and energy through a differential volume of 
fluid, while K(U n ) represents the added contribution 
of the pressure field and centrifugal acceleration to 
the Dalance of radial momentum. The components of the 
vector Uf) are density, axial and radial momentum, 
angular momentum, and total inertial energy. Tne 
symbol S n _i represents tne sum of the contribution 
of tne body force, energy source, and temporal cor- 
relations to the momentum and energy equations, while 
tne subscript n denotes the iteration index of the 
outer loop. The remaining symbol, dv, denotes the 
volume of a differential volume of fluid. Based on 
tne Discussion presentea earlier, the temporal corre- 
lations which are embeaded in 8n_i are simply a 
function of average-passage flow field associated with 
tne neighboring blaae row. Thus, for a single-stage 
configuration, the field equation for the first blade 
row can be written as 

L(uJ ) 1) ) + Jk(U^ 1} ) dv + /s(U n ^) dv » 0 (8) 

while the corresponding equation for the second blade 
row is 


L(uJ. 2) ) + j K(U< 2 >) dv ♦ j S(U n ^) dv = 0 (9) 

Tne superscript (1) and (2) respectively denote 
tne variables associated with the first and second 
blaae row passage flow fields. Next we multiply both 
Eqs. 18) and (9) oy an operator A, which forms the 
axisymmetric average of its argument. This is equiv- 
alent to averaging tne three-dimensional equations of 
motion (i.e., Eqs. (8) and (9)) over the tangential 
direction. For tne first Dlade row, the axisymmetric 
average of the operator L is 

AL(l^ 1) ) » L^'IAUJ 1 )) ♦ j S(U*; 1) ) Ad v (10) 

wnile for the second row 

AL(U^ ) ) - L^IAU^) + /s(U< 2 >) Adv (11) 

In Doth of these expressions, the operator 
denotes the axisymmetric counterpart of L. Tne 
axisymmetric average of the combined integrals which 
appear in Eq. (8) is 

a[/k(uJ 1) ) dv ♦ /s(oJ 2 }) dv] . 

«] + fw£\) Adv (12) 

Similarly, tne axisymmetric average of the combined 
integrals in Eq. (9) is 

a[/k( u J 2) ) dv ♦ /s(U^|)dv] = 

K(U^ 2) ) dv ] + /s(U^j) Adv (13) 



These last results follow because S is independent 
of tangential position. Based on the above equations, 
the axisymmetric average of Eqs. (8) and (9) may be 


expressed as 
l (AX) ( AU U)j + j 

(a[k(uJ 1) ) dv] 



•J 

[s(uj x) ) ♦ s(uJ 2 J 

) ] Adv - 0 

(14) 

L (AX) ( Au (2)) + J 

[a {k(u|; 2) ) dv] 



•/ 

[mu; 2) ) * S(U^J 

)] Adv - 0 

(15) 

Upon convergence of the outer 

loop, Eqs. 

(14) and 


(lb) yield identical solutions for tne axisymmetric 
flow field. In addition, these equations provide a 
means of updating the variables J (U/j) and J (Ufi) 
witnout evaluating the body forces, energy sources, 
and correlations directly. This becomes apparent as 
soon as one notes tnat tne vectors lift and Up ore 
known naving been evaluated in tne inner iteration loop, 

wni le tne quantities /( U S-l> and f ^ U n-1^ are known 

from the previous outer iteration loop. 

This simple strategy for incorporating the 
closure model into a numerical simulation has been 
implemented into the computer code outlined in Ref. 2. 
Tnat code has been used successfully to simulate the 
flow about high-speed counter-rotating propellers as 
illustrated by the results presented in Ref. 2. We 
shall present additional results from that simulation 
as well as that for a high-speed fan stage. 

RESULTS 

Tne model proposed for closing the inviscid form 
of the average-passage equation system was based on 
tne assumption tnat, witnin the confined region of a 
Dlade row, the correlations associated with the blade 
row interaction velocity field are small relative to 
those associated witn the steady aerodynamic blade 
loading. Tne justification for this assumption can 
oe Dased on the argument that the unsteady airload, 
wnich is an indication of the magnitude of the veloc- 
ity component associated with blade row interaction, 
is generally smaller than its time-averaged counter- 
part which serves as a measure of tne magnitude of 
the nonaxisyinmetric component of the average-passage 
velocity field. Data presented in a recent publica- 
tion by Dring et al, (7), shows this to be the case 
in tne midspan region of a turbine stage. Outside of 
tne confines of a blade row, the magnitude of both of 
these velocity fields should be comparable; however, 
their magnitude is small compared to the magnitude of 
the axisymmetric velocity field. As an illustration 
tnat hardware does exist in which one may find regions 
in whicn such flows exist, we present the circulation 
as a function of radius at a number of axial locations 
generated by a high-speed counter-rotating propeller. 
These results are for a flight Macn number of 0.72 and 
an advance ratio for both propellers of 2.8. The cir- 
culation is defined as the integral over a blade pitch 
of tne product of nondimensional radius and nondimen- 
sional tangential velocity. The tangential velocity 
is nondimensional ized by the far-field speed of sound, 
while the radius is rendered nondimensional by the tip 
diameter of tne first propeller. The results pre- 
sented in Figs. 4(a) to (d), are for an axial location 


i 




slightly forward of the first propeller, aft of the 
trailing edge of the first propeller, slightly forward 
i of the second propeller, and aft of the trailing edge 

1 of the second propeller respectively. For each axial 

; location, two plots are drawn. Tne first (solid line) 

| represents the axisymmetric flow field obtained from 

! tne average-passage simulation of tne first propeller. 

The second (dashed line) corresponds to the axisym- 
metric flow field obtained from the average-passage 
simulation of the second propeller, it is quite 
> apparent that both results agree with each other to 

I witnin plotting accuracy and hence are more than 

adequate for assessing blaoe row performance. 

Upstream of the first propeller the circulation must 
i be zero since there is no swirl present in the 

[ incoming flow. At the trailing edge of the first 

propeller, tne circulation is nearly constant over 
the inboard portion, decreasing in a smooth monotonic 
fashion towards zero as the tip is approached (i.e., 
r « 0.5). Thus the aerodynamic loading of the inboard 
region is nearly independent of radius, which implies 
a near-free vortex design. The flow between the two 
propellers would therefore be nearly irrotational . 

The reduction in circulation with radius in the out- 
board region produces a weak tip vortex which convects 
downstream. In the axisymmetric flow representation, 
this tip vortex is smeared into a ring vortex. At the 
leading edge of the second propeller, the circulation 
distribution is seen to be nearly identical to the 
distribution at the trailing edge of the first pro- 
peller. This result further substantiates that the 
flow field between the two propellers is nearly 
irrotational, for in an irrotational unsteady flow 
the time-averaged circulation (or angular momentum) is 
conserved along the axisymmetric stream lines. The 
slight redistribution of circulation that one observes 
in the outboard region is attributed to spanwise mix- 
f ing of angular momentum due to tne tip vortex. 

! Figure 4(d) snows tne distribution of circulation at 

the trailing edge of the second propeller. The second 
i propeller appears to taxe out almost all of tne swirl 

i produced by the first propeller. The change in the 

swirl distribution across the second propeller implies 
that the spanwise aerodynamic loading is nearly uni- 
form over the inboard region of the second propeller. 
The inboard region is behaving as a free-vortex 
design. The results shown in Fig. 4 strongly suggest 
that the closure model developed in this work should 
be applicable to this and similar high-speed counter- 
rotating propellers. This is confirmed by the com- 
parison between the measured and predicted nacelle 
pressure distribution presented in Ref. 2. Further 
comparisons are planned as shown as experimental data 
becomes available. 

An attempt was also made to predict the average- 
passage flow fields generated by a hign-speed fan 
stage. Tne stage chosen was the first of a two-stage 
macnine designed and tested at NASA Lewis (8). The 
computation was performed for an operating point near 
maximum efficiency of tne first stage. This point was 
cnosen to minimize the effect of viscosity on the 
measured flow variables. Tne rotor's rotational speed 
was 80 percent of design and the stage pressure ratio 
and adiabatic efficiency were 1.352 and 0.891 respec- 
tively. The inferred velocity field between the blade 
rows resembled that induced by a free-vortex design 
in the midspan region. The inlet boundary conditions 
in the computation were chosen to produce an inlet 
absolute Mach number and flow angle distribution which 
approximated the measured distributions. At the down- 
stream boundary the nondimensional pressure at the 
hub at the exit of the stator was set equal to the 
measured value, lhe absolute Mach number distribution 


at the inlet to the rotor is shown in Fig. 5, while 
tne relative Mach number distribution across the rotor 
is snowi in Fig. b. The measured results at the inlet 
to tne -otor agree very well witn the predicted 
results, as they snoulo, due to the choice of inlet 
conditions. Tne predicted exit relative Mach number 
distribution appears to be in good agreement with the 
measured results, especially in the midspan region. 

It should also be noted that the predicted relative 
Mach number is less than measured over most of tne 
rotor span. This result is to be expected since the 
blade boundary layers restrict the flow area, thus 
reducing the diffusion capabilities of the rotor. The 
measured distribution also shows the existence of an 
end-wall casing boundary layer which obviously cannot 
be predicted oy the present inviscid analysis. 

The relative flow angle distribution at the lead- 
ing and trailing edge of the rotor was also computed 
ano is shown in Fig. 7 along with the measured dis- 
tribution. This angle is defined as the angle between 
the relative circumferential velocity component and 
the meriodional component. At both stations the pre- 
dicted results appear to be in reasonable agreement 
witn the measurements inboard of the tip region. The 
discrepancy in the tip end-wall region is caused by 
the inability of the present inviscid analysis to 
properly simulate the three-dimensional end-wall flow. 
O'er tne region inboard of the tip, the neglect of the 
influence of viscosity on the simulated axisyinmetric 
flow field produces more turning of the flow than 
experimentally measured, by introducing the effect of 
viscosity into tne current average-passage model 
(whicn includes the outlined closure model), the 
agreement between prediction and experiment should 
improve in tne midspan region. 

Tne next series of results are for the stator. 

Tne absolute Mach number entering and leaving the 
stator is plotted as a function of blade span in 
Fig. 8. The corresponding plots for the aosolute 
flow angle is shown in Fig. 9. This angle is defined 
as the angle Detween the absolute tangential velocity 
component and the meridional component. The agreement 
between the predicted results and measurements appears 
to nave deteriorated from that for the rotor. This 
illustrates the difficulty in predicting multiblade 
row flows. A small error in predicting the perform- 
ance of the first blade row can escalate very quickly 
into a large error in predicted performance of later 
Dlade rows. This problem becomes particularly acute 
whenever there are appreciable regions of flow separa- 
tion in the end-wall region, as appears to be the case 
in tne stator hub region. The poor agreement in the 
stator tip region is attributed to the end-wall wall 
flow induced by tne rotor. To analyze these flow 
regions requires a model whicn incorporates the proper 
ena-wall flow physics. A step in this direction might 
oe maoe Dy including tne effects of viscosity into the 
current average-passage flow solver. An additional 
issue is the development of a closure model for the 
average-passage model applicable to highly rotational 
flows. Kesearcn in botn of these areas is currently 
underway. 

CONCLUSION 

A model was formulated to close the inviscid form 
of the system of equations governing the average- 
passage flow fields for a stage. This model was 
developed so as to insure consistency between the 
average-passage equation system and the axisymmetric 
flow equations. This closure model was used success- 
fully to simulate the average-passage flow fields 


associated with a nigh-speed counter-rotating pro- 
peller. The model was also used in a simulation of a 
nign-speeo fan stage operating near measured peak 
efficiency. This simulation showed tne rotor results 
to be in reasonable agreement with measurements out- 
side of the end-wall region. For the stator the simu- 
lation yielded results wnich were only qualitatively 
correct. The lack of quantitative agreement was 
attributed to neglect of viscosity and the question- 
able applicability of the present closure model to 
the end-wall regions where the flow is known to be 
nignly rotational. Research directed at overcoming 
tnese shortcomings is currently underway. 
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Figure 1. - Flow models for use in turbo machinery design /analysis. 
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Figure 2. - Two stage flow models. 
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Figure 3. - Solution strategy. 
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(c) 2nd Propeller leadinq edge. (d) 2nd Propeller trailing edge. 
Figure 4. - Circulation distribution. 



Figure 5. - Absolute Mach distribu- 
tion (Rotor inlet). 
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Figure 7. - Relative flow angle distribution (Rotor). 
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Figure 8. - Absolute Mach number distribution 
(Stator). 
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Figure 9. - Absolute flow angle distribution (Stator). 
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